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Abstract 

We present here a study of event generators for the hadronic production of Wj . 
We compare a NLO Wj matrix element generator with the leading order calculation 
p<| in Pythia . A matching scheme between a next-to-leading order W^f matrix element 

generator by Baur et. al. and the Pythia parton shower is presented. The NLO 
package produces +0 jet and +ljet final states in the hard scattering and the 
objective is to consistently match these to the initial state radiation from Pythia par- 
ton shower. The proposed methodology preserves both the rate of the hard scattering 
process as well as various kinematic distributions of experimental interest. 
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1 Introduction 



Even though the Electroweak Theory of the Standard Model (SM) has been tested to re- 
markable precision in previous and current accelerator-based experiments, the non-Abelian 
nature of the Electroweak gauge bosons remain to be tested to greater precision. This 
requires a study of the couplings among the vector bosons, the Triple Gauge Couplings and 
the Quartic Gauge Couplings at higher energy scales and the LHC provides an opportunity 
to study the bosonic self-interactions. The cross-section of some of the diboson produc- 
tions, e.g. that of is large enough for a study of pp — > to be feasible at the low- 
energy early running of the LHC with data corresponding to ~100 pb _1 . In preparation for 
analysis of data from the LHC, we present here a study of event generators for simulating 
pp — > event production at the LHC. 

There are broadly two classes of event generators: one which uses exact matrix element 
calculation up to a fixed order and are used for processes involving a few high particles 
(matrix element or ME generators) and the other which uses an approximate scheme to 
generate events up to all orders and involving many particles, most of which are relatively 
soft and collinear (parton shower or PS generators). Event generation also involves a 
hadronization procedure to combine partons from hard scattering and underlying events 
into hadrons which are observable in the experiments. For simulating the hard interaction 
which produces a W and a 7 in proton-proton collision, a matrix element calculation is 
preferred over a parton shower generator. Nevertheless, the parton shower is required 
for providing the countless softer emissions that accompany a hard scattering as well as 
simulating the underlying events and finally hadronizing the product quarks and gluons 
into colour-singlet hadrons, thereby providing a complete picture of an actual event. 

The most popular event generators until now have been of leading order accuracy. How- 
ever, at the LHC, we believe that higher order corrections will become important and may 
lead to event topologies very different from those seen at previous hadron colliders. In view 
of this, we explore the possibility of using a next-to-leading order (NLO) event generator 
for the W"f production process. We choose Baur Wgamma_Nlo [lj as a dedicated NLO 
matrix element generator and compare it to Pythia [2] which is a common general purpose 
leading order generator. Pythia also provides a showering and hadronization mechanism 
for producing the full event. 

The photon transverse momentum spectrum (pj, ) is the most important observable 
parameter in the study of production as this variable is most sensitive to the nature 
of the WW7 vertex coupling. We make a detailed study of p\ from the various source: 
qq' — > Wj hard scattering; photons from initial state radiation from the incoming quarks 
(ISR); photons from the final state radiation from charged leptons from W-decay (FSR); 
and photons from an event with multiple hard scattering, e.g. with qq' — > W"f and qq — > g-y 
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occurring in the same proton-proton collision. 

Eventually we propose a scheme for combining events from the NLO calculation of 
Baur Wgamma_Nlo with the parton shower (which is at leading order) from Pythia 
which takes into account the "double counting problem" of matching a NLO calculation 
with a leading order shower. 

2 W'y production at hadron colliders 

In hadron colliders, the Born level Feynman graphs for Wj events, with W decaying to 
leptons, are shown in Fig. [I] Figure [2] (left) shown the production of W-boson which 
subsequently decays to a lepton and a neutrino with the lepton emitting a photon. The 
quark-gluon fusion process, shown in Fig. [2] (right), where an incoming quark and gluon 
produces a W-boson and a quark in the final state with the quark radiating a photon is 
also present. The FSR diagram leads to the same final state as a W^y production with 
W-decaying through leptons. Other leading and higher order QCD diagrams are present 
as well, some of which are shown in Fig. [3j and their contribution can be quite dominant 
at W'y productions at high energies. The WW'y vertex appears explicitly in the s-channel 
diagram only and is discussed below. 



Figure 1: Born level subprocesses for W'y production in hadron-hadron collision : s-channel 
process (left), t-channel process (middle) and u-channel process (right). 



Figure 2: Final state radiation from the charged lepton from W-decay (left) and Quark- 
gluon fusion diagram producing a W'j in the final state(right). 
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Figure 3: Higher order QCD diagrams for production in hadron-hadron collision. 
2.1 WW^( vertex and cross-section 

The most general Lorentz and electromagnetic gauge-invariant CP-conserving Lagrangian 
for the WW^i vertex can be written as follows: 
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where A^ and W 1 are the photon and W _ -fields respectively, = d^W u — d v W^ and 
Fp,v = d^A u — d u A^. The variables k and A are related to the magnetic dipole moment nw 
and the electric quadrupole moment Qw of the jy-boson: 
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The vertex term in Eqn. [T] contains both the SM and the non-Standard Model (NSM) 
contributions, the latter being expressible in terms of An = k — 1 and A, both of which are 
zero in the SM. 

The differential cross-section for q^q 2 — > W - '-/, in the Standard Model [3J, with k — 1 
and A = 0, is given by 
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where s, t and it are the Mandelstam variables, gi2 = cosOc for q\q 2 = du and sc and 
t?i2 = sinOc for gig 2 = sm and dc. Q±e is the charge of the quark q\ and Q 2 = Qi + 1. The 
differential cross-section can also be expressed as 



(3) 



where #* is the angle between the W and ci-quark or equivalently between the 7 and 
w-quark in the centre of mass frame of the system. 
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A special feature of the gauge theory is manifested in this distribution due to the value 
of k = 1 and A = in SM. Notably, the differential cross-section da(qiq 2 Wj)/dcos9* 
vanishes at a particular angle 9*: cosO* = —1/3. This phenomenon, the radiation amplitude 
zero (RAZ), is possible only for SM values of k and A. This feature can be attributed to the 
expression [Qi + 1/(1 + t/u)} 2 which vanishes for t*/u* = —(1 + 1/Qi), or in other words, 
for cosO* = —(1 + 2Qi) corresponding to the charge of the quark Q\ = —1/3. 

The zero occurs due to destructive interference of radiation patterns in gauge theory 
tree level amplitudes for the emission of massless gauge bosons [I], [5], [5]. Any anomalous 
moment resulting in different values of the coupling, would destroy the occurrence of the 
zero. Consequently the precise measurement of the radiation amplitude zero would serve to 
establish the Standard Model and also to search for Beyond Standard Model (BSM) physics. 
The anomalous couplings can be tested in such a situation by measuring the production rate 
at high value of photon transverse momentum. The DO experiment has recently reported 
a study of the radiation amplitude decay at the Tevatron [7j. In the present study, we deal 
only with the Standard Model WW7 vertex parameters, i.e. Ak = A = 0. 

2.2 Event generators 

In preparation for data analysis for study of Wj events produced at the LHC, we would 
like to choose an event generator that has all the features of Wj production in proton- 
proton collision at the energy of the LHC. W'y production. The common workhorse in the 
CMS collaboration for generating events is the Pythia event generator. It is a general 
purpose generator containing matrix element calculation of many 2—7-2 and 2— >1 processes 
as well as showering and hadronization scheme. The calculation for W^f process is present 
up to Born level, i.e. it includes all graphs shown in Fig. [TJ QCD radiation can be added 
by parton showers as well as photon bremsstrahlung from the outgoing charged lepton 
from W-decay using parton shower FSR. Yet, given that the probability of a hard parton 
emission accompanying the W'-f production at the LHC and the possibility of the WW^ 
vertex containing anomalous coupling parameters, we would like to choose an NLO matrix 
element generator with anomalous couplings. We have chosen the Baur Wgamma_Nlo 
generator and will present here a comparison of Baur Wgamma_Nlo with Pythia . 

3 Comparison of Baur Wgamma NLO and PYTHIA 6 

3.1 Baur Wgamma NLO at the Tevatron and LHC energies 

We first put forward a comparison of Baur Wgamma_Nlo at the energy of the Tevatron 
(1.8 TeV) and at the LHC (10 TeV). The cuts on various Baur Wgamma_Nlo parameters 
for the Tevatron is listed in Table [T] and those for the LHC are in Table [2] 
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Table 1: Generator level cuts for Baur Wgamma_Nlo for generating events at 1.8 TeV 
proton-antiproton collisions. 
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Jet pt 


1.0 GeV 


Photon rapidity 


1.0 


Charged lepton rapidity 


2.5 


Jet rapidity 


10.0 


AR(j, lepton) 


0.7 


Cluster(W,7) transverse mass 


90 GeV 


Soft divergence parameter 


0.01 


Collinear divergence parameter 


0.001 



Fraction of hadronic energy in a cone around the photon 0.15 



Table 2: Generator level cuts for Baur Wgamma_Nlo as used for generating LHC events 
at lOTeV. 



Parameter 


Cut 


Photon pt 


5 GeV 


Charged lepton p T 


5 GeV 


Neutrino pt 


5 GeV 


Jet pt 


5 GeV 


Photon rapidity 


8.0 


Charged lepton rapidity 


8.0 


Jet rapidity 


8.0 


AR(j, lepton) 


0.05 


Cluster(W,7) transverse mass 


10 GeV 


Soft divergence parameter 


0.01 


Collinear divergence parameter 


0.001 



Fraction of hadronic energy in a cone around the photon 0.15 



Fig. [J] (left) shows the p T spectrum from BAUR WGAMMA_NLO at the Tevatron. The 
two curves shown are for the Born level cross-section and the next-to-leading order cross- 
section. The k-factor is shown in Fig. [i] (right). 

At the energy of the LHC, we expect an overall enhancement of the cross-section, as 
compared to the Tevatron energy, as shown in Fig. [5] (left). Further, the k-factor value 
(Fig. [5] (right)) is also larger and is no longer a constant but increases with increasing value 
of the photon transverse momentum. Thus, for high p T we see that the NLO corrections 
become very important. 
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Figure 4: Photon p T from Baur Wgamma.Nlo at ^Js = 1.8TeV (left) and the k- factor, 
as defined in the text (right). 



> 




Figure 5: Photon px from Baur WgammaJMlo at LHC centre of mass energy, lOTeV 
(left) and the k-factor at that energy (right). 

3.2 PYTHIA 6 matrix element vs Baur Wgamma_NLO Born 
level 

We compare Baur Wgamma_Nlo with Pythia 6 at 10 TeV proton-proton collision. 
Using Pythia , we have the following options to produce W^y events : 

1. As mentioned before, Pythia can produce events with full matrix element cal- 
culation of the Born level Feynman diagrams. We compare this Pythia 2— >2 process 
withe the events obtained from using Baur Wgamma_Nlo 's calculation at Born 
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level. The only difference between Baur Wgamma_Nlo ad Pythia at this stage is 
that Baur Wgamma_Nlo includes the decay of the W-boson to lepton in its matrix 
element calculation but in Pythia the W decay is treated separately. 

2. Additionally, we can generate inclusive W events in Pythia and add photon emissions 
using the PYTHIA parton showers. It is worthwhile to compare the photon transverse 
momentum spectrum of photons generated from the parton showers to that from 
Pythia 's matrix element calculation and also to Baur Wgamma_Nlo 's. 

The Baur WgammaJNlo parameters from Table [2] are used. For Pythia 6 matrix 
element (option 1 above), we switch on the Wj production process. ISR and FSR were 
switched of and so was the primordial kj- of the partons inside the protons. 

Figure [6] (left) shows the jP T spectrum from Baur Wgamma_Nlo Born level and 
Pythia matrix element. We see a good agreement between Pythia and Baur Wgamma_Nlo 
at this level. Figure [6] (right) shows the pj~ spectrum from Baur Wgamma_Nlo NLO 
calculation and compares it to Pythia ISR photons from parton showers in inclusive W- 
boson events. The ]P T distribution from Pythia matrix element is also plotted on the same 
canvas for perspective. We see that NLO effects do play a major role in altering the shape 
of the photon transverse momentum distribution. However, the parton shower-generated 
photons mimic the spectrum from the NLO events quite well. This is because the parton 
showers are supposed to approximate the calculation of a process up to all orders. The only 
problem with the parton showers is that they are mainly designed for soft and collinear 
emissions and so events at the high end of the transverse momentum distribution are very 
infrequently generated. 

Also parton showers in Pythia 6 do not have anomalous couplings and the radiation 



zero. As mentioned in Section 4.1, Pythia 8 parton shower also contains the radiation 
zero feature, included as matrix element correction to the QED component of the parton 
shower. 

3.3 PYTHIA 6 effective "k-factor" 

Just as the Baur Wgamma_Nlo generator has leading order QCD diagrams with a quark 
or gluon emission contribution to a boost to the Wj system, the Pythia ISR QCD radi- 
ation also serves the same purpose and contains showers up to all orders although to some 
approximation. This leads to defining an effective k-factor for the Pythia ISR as well 
which is given by the ration of Pythia cross-section for the W^f process with ISR turned 
on to that with ISR turned off: 

da ISR ON 

Pythia "k - factor" = ^J^uff 
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Figure 6: Photon transverse momentum spectrum from BAUR WGAMMA_NLO Born level 
and Pythia ME (left) and the spectrum from Pythia matrix element calculation 
without PS compared to Baur Wgamma_Nlo NLO calculation (right). Superimposed 
on the same canvas (right) is the photon p t spectrum from Pythia ISR in an inclusive 
W-boson production. 



Figure [7] (left) shows the plot of this quantity. The ISR QCD emission, as anticipated is 
much softer than the matrix element calculation of Baur Wgamma_Nlo for a single jet 
emission. Figure [7] (right) shows the transverse momentum of the W'y system from Baur 
Wgamma_Nlo NLO and Pythia ISR. These two plots show that the QCD radiation 
from Pythia 6 parton shower is comparatively softer than the hard emission from Baur 
Wgamma_Nlo matrix element calculation; the Baur Wgamma_Nlo k-factor, as shown 
in Fig. [5] (right) has a higher value compared to the Pythia 6 equivalent. 

4 PYTHIA 6 vs. PYTHIA 8 and an improved parton 
shower description 

With Pythia 8 [8] , we have a a few advantages over Pythia 6. Pythia 8 gives an 
improved description of the parton shower description, multiple hard-scattering and a code 
written in C++. 

The parton shower of Pythia 8 is compared to that of Pythia 6 in Fig. [8j Figure [8] 
(right) shows the parton shower photon pt spectrum from PYTHIA 6 and PYTHIA 8 using 
their default settings. The agreement is not perfect and a complete agreement between both 
the Pythia versions are achieved by setting the Pythia 8 parton shower setting to the 
default vale of Pythia 6 as shown in Fig. [9j This comparison is important because 
in future we use Pythia 8 for our matching parton showers to Baur Wgamma_Nlo 
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P^GeV) pw+pj, (GeV) 

Figure 7: The ratio of spectrum from Pythia matrix element with and without ISR. 
This gives an effective "k-factor" for the Pythia parton shower which is supposed to take 
into account higher order emission diagrams for the process (left); and the plot of the 
transverse momentum of the Wj system from Baur Wgamma_Nlo NLO and Pythia 
inclusive W events with ISR. This gives the amount of boost received by the W7 from QCD 
radiation (right). 



-generated events and a consistency check is required. 




Figure 8: Comparison of matrix element photon from PYTHIA 6 and PYTHIA 8 (left) 
and comparison of parton shower photons from PYTHIA 6 and PYTHIA 8 (right). 

4.1 ISR, FSR and Multiple Scattering 

Pythia 8 includes an improved description of the parton shower which agrees more closely 
with the NLO calculation from Baur Wgamma_Nlo . The parton shower photon in 
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Figure 9: Comparison of matrix element photon from PYTHIA 6 and PYTHIA 8 (left) 
and comparison of parton shower photons from PYTHIA 6 and PYTHIA 8 (right) with 
same settings for both Pythia 6 and for Pythia 8. We see that the agreement is better 
compared to that in Fig. |8j 

Pythia 8 now also contains the radiation zero which is present in the qq — > W'y calculation. 

Sources of photons, as mentioned before, could also be the FSR photon from the outgoing 
charged lepton if the W decays via leptons, as in Fig. [2] (left). At the detector level, this 
yields the same final state as processes of Fig. [T] and so the events containing FSR need 
to be eliminated to select the actual events. Baur Wgamma_Nlo does not contain this 
diagram Fig. [2] (left) and Pythia adds this diagram through the parton shower evolution 
of the FSR. Thus using Pythia it is necessary to understand the contribution of these 
photons to the overall photon distributions. 

Similarly, a second hard process, e.g. qg — > jq may occur with a W-boson production 
via qq' — > W and will have all the characteristics of a Wj event. The new Pythia 8 
provides the setting to study this contribution as well. 

Figure [10] shows the photon spectrum from both the above sources and compares it 
with the photons from PYTHIA 8 ISR and matrix element from BAUR WGAMMA_NLO 
. The FSR photons do not survive above 40 GeV which is the kinematic cutoff and a 
reasonable cut of the pj. on events from the detector should remove almost all of these. 
Interestingly, the contribution from multiple hard scattering is also insignificant compared 
to the photon counts from the hard scattering process and ISR and this too is confined to 
a low region of p^ . 
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Figure 10: Comparing photon from PYTHIA W'y matrix element, from ISR shower and 
both ISR and a second hard interaction producing a photon in the same event. 



5 Matching 

This section is based entirely on the C++-based version Pythia 8. Consequently, any 
reference to Pythia would mean Pythia 8.120, unless otherwise stated. 

A collision of two protons involve a hard scattering part which produces the signal of 
interest as well as many other soft collisions for the other constituents of the proton, the 
quarks and gluons (collectively called partons), which need to be taken into account to make 
a complete description of the event. The debris of the collisions, both hard and soft, are 
never actually seen in the detector. Instead, one sees other particles which are composites 
of these fundamental partons which occur through he process of hadronization. This is a 
unique feature of the interaction that governs particles like protons and its fundamental 
constituents the partons. The evolution of the partons into more composite particles, known 
as hadrons is called hadronization and what we can detect are these hadrons. It is obvious 
that for claiming to model an event as a real life scenario, a Monte Carlo program needs to 
take care of the hadronization of quarks and gluons. Baur being a dedicated event generator 
for W+photons events does not take into account either the soft part of a collision nor of 
the hadronization. On the other hand, Pythia being a general purpose event generator, 
provides a full description of the event with the hard scattering, underlying events, parton 
showers adding the ISR and FSR and finally the hadronization of the products. 

We would like Baur Wgamma_Nlo to simulate the hard scattering two partons to 
produce a Wj event with W decaying to, say a muon and muon-type neutrino. After this, 
we will use Pythia to take care of ISR and FSR through the parton showers and then also 
simulate the underlying events and the hadronization of the final products. However the 
following problems are encountered if we naively let Pythia parton shower act upon the 
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events from BAUR WGAMMA.NLO : 

Since Baur Wgamma_Nlo is a matrix element event generator, it produces each event 
associated with a certain weight: i.e. it samples the phase space available uniformly and 
assigns a weight to each event corresponding to the probability of that region of phase 
space getting populated. This poses a problem because we want events to be generated 
with unit weights. Second, there is the problem of "double counting" in certain regions of 
phase space, which we discuss below: 

BAUR WGAMMA_NLO produces two kinds of events: events with no partons in the final 
states, with only the W and the photon with W — >■ /m> and events with a parton in the 
final state apart from the W decay products and the photon. The former will be called 

3- body events and the latter, 4-body events. These 4-body events are the exclusive 
+ljet events, given some suitable definition of a jet. The Baur Wgamma_Nlo program 
calculates the cross-section of both the 3-body (0-jets) events as well as 4-body (1-jet) 
events. Addition of parton showers by Pythia should preserve these cross-sections to first 
order. But it may so happen that during the showering of some of the 3-body events, the 
parton shower generated by Pythia agrees with the definition of the jets for the 4-body 
events from Baur Wgamma_Nlo . Thus these 3-body events would be considered as 

4- body (1-jet) events after the showering, thereby altering the exclusive cross-section of 
both W'-f and W'-f +ljet events. 

We have demonstrated a proof-of-principle matching scheme for Baur and Pythia and 
further details like looking that the parton spectrum and cuts on the parton (jet definition) 
has to be optimized. 

5.1 The matching algorithm 

1. Baur WGAMMAJNLO produces 3-body final states (//, v and 7) and 4-body final 
state (/!, v, 7 and jet (q or g)). 

2. The 4-body events lack a Sudakov form factor (the probability of no emission which 
take into account the effect of virtual loops) which requires the following: 

• Project onto a 3-body, by assuming that the outgoing parton can be emitted from 
either incoming partons (flavours permitting), with relative weights by splitting 
kernels and parton densities. So we assume that the 4-body state never had a 
parton emitted and all the kinematics are recalculated based on that. This gives 
the projected 3-body event. 

• Shower the resultant 3-body, compare p^ ower at the first ISR branching with the 

parton • • • 1,11 1 

m original 4-body. 

• If p^ ower > p^ rton then the event is reclassified to 3-body; move to step 3 below. 
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• If p^ ower < p T arton then the original 4-body event is now showered and p^ ower is 

j 'ii parton 

compared with p T 

• If p^ tower > p^ Lrton then go back one step. 

• Continue with rest of shower to give a complete event. 

3. Shower the 3-body events; compare p^ ower with p^P arate after first ISR branching. 
^separate can ^ g cons j ( j ere( j £ De boundary between the ME calculation's regime 
and that of the parton shower's. 

4. If p^ lower > pf^P arate then stop any further shower evolution and go back one step. 

5. Continue with the rest of the shower. 

5.2 Matching results 

According to the above-mentioned algorithm, we first project 4-body states on to 3-body 
states and try to determine the Sudakov form factor for the QCD emission which is originally 
from the matrix element calculation. 

Figure [Tl~| (left) shows the transverse momentum of the jets after the Sudakov correction 
and the px of the jets matched after branching. The delta-function at the origin correspond 
to those events which are reclassified as 3-body events after the Sudakov correction. 

As a second step to the parton matching scheme, the events are showered and the 



first ISR pt is plotted in fig 11 (right). Since according to our algorithm, we are clearly 
demarcating a px region above which matrix element calculation is valid and below which 
is the parton shower regime (at 5 GeV for the present study) we see that the first ISR 
emission from 3-body events are always confined to below 5 GeV. But this is not so in 
the case where the parton showering is only required to be softer than the matrix element 
parton. Hence for 4-body events, we get a long tail for events with highly energetic matrix 
element parton. 



Figure |X2| shows the shows the distribution of the transverse momenta of the system 
as produced by Baur Wgamma_Nlo and after the Pythia showering. The red histogram, 
depicting the distribution from Baur Wgamma_Nlo shows many events with p^ 1 equal 
to zero which are the 3-body events. The non-zero values correspond to events with a parton 
in the final state. In the green histogram, we see that the kink after the zeroth bin fills up 
due to the boost from Pythia ISR. The area under both these curves however remain the 
same indicating that the exclusive cross-section of the 1-jet events remain conserved after 
the parton shower. 

While being suitable for a simple event topology with one jet, the matching scheme does 
not require any modification of the matrix element calculation of Baur for the Sudakov form 
factor. 
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Figure 11: Transverse momentum of the partons after Sudakov correction and matching 
(left). Events which were originally of the type 4-body and subsequently gets reclassified as 
3-body after the Sudakov correction end up with no partons in the final state. These events 
show up as the delta-function peak in the Oth bin of the red histogram. The Transverse 
momentum of the partons from 3 and 4-body events after 1st ISR emission are shown on 
the right. 




Figure 12: Transverse momentum of the W'y system before Pythia showering and after 
Pythia showering. 
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6 Summary 



The preparation for data analysis of high energy physics experiments involve the use of 
Monte Carla techniques to simulate generation of events one wishes to study and also simu- 
lation of detector equipment and interactions of the produced particles in the detector. As 
a preparation for studying the production of W7 events at LHC proton-on-proton collision 
using the CMS detector, we wish to accurately model the production mechanism of W+7 
at the LHC. In this note, we have made a comparison of a well-established general purpose 
event generator, Pythia with the dedicated Baur WgammaJNlo generator. Pythia 
has been found wanting on several issues regarding W^y production mechanism and Baur 
Wgamma_Nlo too has been found to lack the total description of a realistic hadron col- 
lision scenario leading to W'y production. The most effective strategy, thus would be to 
combine the two event generators to the best of our advantage and in the later section 
of this note, we have outlined a home-brewed approach to this combination, the so called 
"Matching strategy". Given the individual performance of Baur and Pythia, a successful 
implementation of the matching strategy will pave the way for event generation taking into 
account aspects of all physics involved, at various scales, in hadronic production of W-y at 
the LHC. 

References 

[1] U. Baur, T. Han, and J. Ohnemus, Phys. Rev. D48 (1993) 5140-5161, 
[hep-ph/9305314| . 

[2] T. Sjostrand, S. Mrenna, and P. Skands, JEEP 05 (2006) 026, |hep-ph/0603175] . 
[3] R. W. Brown, D. Sahdev, and K. O. Mikaelian, Phys. Rev. D20 (1979) 1164. 
[4] R. W. Brown, K. L. Kowalski, and S. J. Brodsky, Phys. Rev. D28 (1983) 624. 
[5] S. J. Brodsky and R. W. Brown, Phys. Rev. Lett. 49 (1982) 966. 
[6] M. A. Samuel, Phys. Rev. D27 (1983) 2724-2731. 

[7] V. M. Abazov et. aL„ DO Collaboration Phys. Rev. Lett. 100 (2008) 241805, 
|0803.0030] . 

[8] T. Sjostrand, S. Mrenna, and P. Z. Skands, Comput. Phys. Commun. 178 (2008) 
852-867, |0710.3820]. 



16 



